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Abstract. 

Two-level atoms interacting with a one mode cavity field at zero temperature 
have order parameters which reflect the presence of a quantum phase transition at 
a critical value of the atom-cavity coupling strength. Two popular examples are the 
number of photons inside the cavity and the number of excited atoms. Coherent 
states provide a mean field description, which becomes exact in the thermodynamic 
limit. Employing symmetry adapted (SA) SU(2) coherent states (SACS) the critical 
behavior can be described for a finite number of atoms. A variation after projection 
treatment, involving a numerical minimization of the SA energy surface, associates 
the finite number phase transition with a discontinuity in the order parameters, which 
originates from a competition between two local minima in the SA energy surface. 



PACS numbers: 42.50.Ct, 03.65.Fd, 64.70.Tg 
1. Introduction 

Cavity quantum electrodynamics describes a system containing elements which can 
absorb and emit radiation interacting with quantized radiation. The simplest example 
is that of one atom interacting with a single cavity mode, which exhibits the presence 
of quantum revivals and squeezing [TJ |2j. For N non-interacting two-level atoms [3], 
above some critical coupling strength, all atoms are predicted to absorb and emit 
light collectively and coherently, in a superradiant mode [4j. This mode has been 
observed in a driven system [5] which, under certain assumptions, is described by the 
same dynamics [6]. An analytical description, based in a truncated Holstein-Primakoff 
expansion and valid in the thermodynamic limit, when the number of atoms goes to 
infinity, has been widely used [SI El El E] , although it is divergent at the critical point. 
Finite size corrections have been studied in Ref. [10]. A similar behavior is found in 
superconducting circuit QED systems [5], where the number of elements is small and it 
is important to take explicitly into account these finite size effects. 

The description of the phase transition at zero temperature and for a large number 
of atoms can be performed analytically through the mean field formalism. The radiation 
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is described with harmonic oscillator, Heisenberg-Weyl (HW) coherent states (CS) 
[T2| [13]. The atomic sector can be described with HW-CS, employing the Holstein- 
Primakoff realization of the angular momentum operators [7J [TT], or with spin SU(2) 
CS [14J. In both cases the energy per atom, the number of photons per atom and 
the fraction of excited atoms are very well described in the large TV limit, while other 
observables are misrepresented due to the violation of the Hamiltonian parity symmetry 
in the description of the ground state wave function in the superradiant phase. This 
difficulty can be overcome by employing symmetry adapted coherent states (SACS), 
which by construction restore the parity symmetry [151 EE E] • 

While SACS allows for a better description of the ground and first excited states 
than CS, the energy surface obtained with SACS is more complicated, preventing a direct 
analytical treatment of the critical points. Employing the CS expressions obtained for 
them in the SACS description provides wave functions which have good overlap with 
the exact wave function calculated at finite N (but for a small vicinity around the phase 
transition), allowing to express the expectation values of matter and field observables 
with elegant analytic forms in the superradiant regime [18] HU |20l |2U [22] . Many of 
them coincide with the mean field ones for large N, while some others, noticeably the 
fluctuations, differ. This analytic description fails, however, to describe effects of order 
1/JV in the normal phase at finite N, which are absent in the mean field description. 

To go beyond this point, we have performed a variational treatment of the SACS 
energy surface. In what follows we present the details of this formalism, which exhibits 
a phase transition at finite N at values of the coupling constant which are larger 
than those found in the thermodynamic limit, and which has non-trivial values of the 
order parameters, in agreement with the numerical calculations. These results can be 
applied to the description of the phase transition at large but finite N, where numerical 
simulations are difficult. 

2. The Dicke Model 

The Dicke Hamiltonian allows the description of N two-level atoms, with energy 
separation u^, interacting collectively with a one mode radiation field, whose frequency 
(times h) is employed as the unit of energy along this work. Photon creation and 
annihilation operators are represented by the operators a) and a, the atomic relative 
population operator by J z , and the atomic transition operators by J±. The Hamiltonian 
reads 

H D = a) a + oo A J z + ^= (a* + a) ( J_ + J + ) , (1) 

where 7 is the (adimensional) atom-field coupling parameter. 

The eigenvalues j associated with the amplitude of angular momentum operator J 2 
are fixed as j = N/2 along this work, where we are restricting our analysis to the fully 
symmetric representation of the atomic states, which are enough to describe its collective 
behavior. When there is no interaction the number of excited atoms is j + m, where 
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m is the eigenvalue of J z ; the number of photons inside the cavity is u, the eigenvalue 
of the photon number operator a^a; and the excitation number is A = j + m + u, 
the eigenvalue of the excitation number operator A. A is a conserved quantity in the 
Jaynes- and Tavis-Cummings models [23j[2l], where the rotating wave approximation is 
employed [J71 125], but the counterrotating terms in the Dicke Hamiltonian (a* J + , a J J) 
mix them. This has far reaching consequences: there are no finite subspaces where the 
Dicke Hamiltonian can be exactly diagonalized, it is not integrable, and around the 
phase transition it exhibits the presence of quantum chaos (26J, E7J . 

There is, however, a remnant symmetry: the parity symmetry operator e l7rA 
commutes with the Dicke Hamiltonian, allowing the classification of its eigenstates, 
as they belong to the symmetric or to the antisymmetric representations of the cyclic 

group c 2 nana nu. 



3. The mean field analysis 

As mentioned in the introduction, the mean field treatment allows for the appropriate 
description of many intensive observables in the thermodynamic limit. In this section we 
briefly review the main results, which have already been reported [TS1 [JJJ]- The starting 
point is a variational state built as a direct product Heisenberg-Weyl coherent states 
|a) [121 H3] for the photon sector and SU{2) or spin coherent states \() [THUS] for the 
particle sector, i.e., \a, () = \a) <S> |0- They are 

00 a u 

\a) = (exp exp(aa^)|0) = (exp ^2 — ^= \v), 

z,=ovV! 

10 =(l + IC| 2 )^ /2 exp(CJ + )|j, -j) (2) 

m=—j \ J / 

where the ket |z/) is an eigenstate of the photon number operator, and \ j, m) is a Dicke 
state. This trial state describes = 2j particles distributed in all the possible ways 
between the two levels, and up to an infinite number of photons in the cavity. 

Expressing the internal variable a in terms of the quadratures (q,p) associated 
with harmonic oscillator realization of the field, and ( in terms of the angles (8, 4>) 
which determine its stereographic projection on the Bloch sphere, 

a = ^(q + ip), C = tan exp {-i 0) , (3) 

the expectation value of the Hamiltonian with respect to the trial states is obtained. 
From it, the energy surface is defined dividing by the number of particles. It is 

E(q, p, 9, 4>) = ^ (f + q 2 ) - cos 9 + sin 9 cos <p , (4) 

It depends on two control parameters (u, 7) and four order parameters (q,p,9,(j)). 
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The catastrophe formalism [28] allows the obtention of the critical points 
Qcj Pc C , 4>c and the associated separatrix, the geometric place in the control parameter 
space where a catastrophe, i.e. a sudden change in the ground state of the system, 
reflected in the order parameters, takes place. They are the solutions of the system of 
equations 

dE 

— = 0, u = q,p,9, 0. (5) 

Those which correspond to minima in the energy surface are [TU [19] 

C = O, <? c = 0, Pc = 0, for |7| < 7 C 

9 C = arccos(7 c /7) 2 , q c = -2 yf] 7^1 - (7 c /7) 4 cos0 c , p c = 0, for | | > 7 

The condition I7I < 7 C defines the normal phase, where the ground state has zero 
photons and no excited atoms. The superradiant phase, described in the second row, 
assumes oja > 0. Introducing the variable x = 7/70 the energy at the minima takes the 
values 



(6) 



c 



-^normal 2A^7 C , -^superradiant N / y c x (l + x . (7) 
4. Symmetry- Adapted Coherent States (SACS) 

The variational states which preserve the symmetry of the Dicke Hamiltonian can be 
built from even and odd parity coherent states |a) <8> |0. They have the form [18j [19] 

|«, C)±=A4(|«)®|C)± \-a)®\-()) , (8) 



with normalization M- 



± 



j ^ , = 2 ( i = t «^(- 2 H*)(y?pj j ■ < 9 » 

The above expressions can be used to evaluate the expectation value of the Dicke 
Hamiltonian. When divided by the number of particles, it represents the energy surface 
with respect to the symmetry adapted states (SAES), which has the form [19] 

2 1 



\ (p 2 + <r 



E±= ± - iv z + a z ) { 1 - 



N [ +1 tan 2 ^ cos^ ] 

uja \ (cos 6M ± 



2 [ v ' l±e±(f 2 +r)(cos0)^ J 

/ — — ( ±p tan^ sin0 + qe p2+q2 sin# cos0 (cos^) - ^ 1 
+ V2iY7 \ e^ 2 (cos^±l J" (10) 

In the limit — > 00, E + reduces to Equation Q, when |cos^| 7^ 1. However, it is 
important to stress that the analysis can be carried out for any value of N, and in this 
contribution we work at finite N. 

In our previous works we have used the SACS, Equation Q, with the critical points 
of the standard CS given in (J6j) . It allowed us to obtain analytical expressions for all 
the observables of interest |T9]. In the limit iV — > 00 they coincide with those obtained 
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with the coherent state, with the exception of (Ag) 2 and (AJ X ) 2 . These two exceptions 
arise because (q) = (J x ) = for the SACS but not for the CS. As expected, the SACS 
results are in close agreement with the numerical calculations for finite N [IS] . 

There are, however, clear limitations in this mixed approach, where the coherent 
states have the parity symmetry restored but the critical points are not obtained by 
minimizing the corresponding energy surface. In the normal phase the solutions are 
trivial but not correct for finite N, because effects of order 1/N are not considered at 
the mean field level. This problem extends into the superradiant phase: in a vicinity 
of the phase transition the functions employed do not provide a minima of the energy 
surface given in Equation (10), because they are not solutions of ([5J [19] . 



5. Full variational treatment 



The exact critical points of the SAES (10) are obtained by applying to it the conditions 
given in ([5]), i.e. its first derivatives with respect to the four internal variables should 
cancel. The partial derivative with respect to p implies p = 0. Using it, the partial 
derivative with respect to <fi leads to sin0 = 0. For u>a > 0, the minima occurs for 
= 0. Introducing z = cos 9 the other two equations which must be satisfied, for the 
even parity case, are: 



dE+ 1 _ 2q2 



dq z 



e 



qz 1+4j + e 2q z[q + 2Jj -jz 2 -f 



+2e q z 3 [q z + zJj - jz 2 ^ + 2q zJj - jz 2 ^ -jq[-l + z u A 



dE, 1 



89 z 26 



-2g 



-jz^VT^uja + e 2q2 (2 x /jqz 3 -f + jz 2 ^^ 2 



.2 



Z*U A 
3/2 



+ e q2 ^fjz 2j ( 2gz 3 7 - Ajqz (-1 + z 2 ) 7 + 2j 3/2 (l - z~ J ~c , 



+JjVl -z 2 (2q 2 z + (-1 + z 2 ) 



0. 



We have explored two ways to find the minima numerically: solving the above 
set of coupled non-linear differential equations, and exploring the SAES to find them. 
Employing Mathematica [2H], we found that the most efficient way was the second one. 
In the following figures the symmetry adapted energy surface is plotted for five values of 
the interaction strength 7, at resonance (ua = 1). On the left hand plots, the horizontal 
axis represents q, and the vertical axis 9. Continuous lines connect regions with equal 
energies, and the minima are shown with large dots 

On the right hand side, the SAES is plotted against q, with 9 varying as a 
linear function of q which crosses both minima. Figure [I] displays the results for 
7 = 0.545, 0.550, 0.552, and figure [2] for 7 = 0.555 and 0.560. It is clearly seen that 
for 7 = 0.545 the SAES has a global minimum around q = 1, and a local one around 
q = 2 is very shallow. It becomes deeper for 7 = 0.550, where the two local minima 
are well developed, but still the first one is the more profound. At 7 = 0.552 the two 




Figure 1. The symmetry adapted energy surface showing the minima for 7 = 
0.545,0.550,0.552 plotted as a function of q and (left) and as function of q (right). 
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Figure 2. The symmetry adapted energy surface showing the minima for 7 = 
0.555,0.560 plotted as a function of q and 9 (left) and as function of q (right). 

minima have the same depth. This is the point where the phase transition takes place. 
The two cases shown in figure [2] display the evolution of the absolute minimum around 
q = 2. 

The phase transition in the Dicke model, in resonance, for finite number of atoms 
N = 20, takes place at a critical value 7 C = 0.552, 10% larger than the critical value 0.5 
in the thermodynamic limit. As expected, with increasing N it tends to this limit. The 
existence of the phase transition for finite values of N has also been found as a singularity 
in the fidelity susceptibility obtained through exact diagonalization [30l [31]. In the 
Dicke model they are similar to those found using SACS [32J. There are, however, some 
important differences. On one hand, the SACS approach provides a description of all 
observables and of the phase transition for any number of atoms, avoiding the difficulties 
inherent to the diagonalization. Two order parameters (the number of photons inside the 
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cavity and the fraction of excited atoms) have finite values (g c /2 and cos 6^, respectively) 
in the normal phase, improving in this way the mean field treatment. On the other 
hand, the jump in the critical variables q c , 9 C reflects directly in a jump in these order 
parameters |32j . While it serves as an indicator of the presence of a phase transition at 
finite N, this discontinuity is spurious, in the sense that it is not present in the exact 
numerical calculations. At the same time, there is a universal relationship between the 
two order parameters which is exact, and is fully satisfied by the present description 
[32J. A detailed analysis of the critical strengths 7 C for a finite number of atoms will be 
presented elsewhere [33] . 

6. Conclusions 

The phase transition in the Dicke model, for a finite number of atoms N, has been 
presented employing symmetry adapted coherent states. Analytic expressions for the 
symmetry adapted energy surface were presented. It was minimized numerically. It was 
shown that around the phase transition it has two minima, and the competition between 
them determines the absolute one. It allows for a precise determination of the critical 
strength for any given N, which tends to its mean field value in the thermodynamic 
limit. It provides finite values for the number of photons inside the cavity and the 
fraction of excited atoms in the normal phase, improving the mean field results, but at 
the price of having a discontinuity at the phase transition, which is erased in the limit 
iV — > oo. The Dicke model keeps offering fascinating new experimental and analytical 
results. 

This work was partially supported by CONACyT-Mexico, and DGAPA-UNAM 
under project IN102811. 
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